ACCEPTED BY ApJ LETTERS JUNE 23, 2009 

Preprint typeset using I^T^X style cmulateapj v. 10/09/06 



COLLISIONAL STRIPPING AND DISRUPTION OF SUPER-EARTHS 
Robert A. Marcus 1,a , Sarah T. Stewart 2 , Dimitar Sasselov 1 , Lars Hernquist 

Astronomy Department, Harvard University, Cambridge, MA 02138 and 
2 Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA 02138 
accepted by ApJ Letters June 23, 2009 



OS 

o 
o 

(N 



(N 

w 

on- 
6 



> 

q 

o 

o> 
o 



X 



ABSTRACT 

The final stage of planet formation is dominated by collisions between planetary embryos. The 
dynamics of this stage determine the orbital configuration and the mass and composition of planets in 
the system. In the solar system, late giant impacts have been proposed for Mercury, Earth, Mars, and 
Pluto. In the case of Mercury, this giant impact may have significantly altered the bulk composition 
of the planet. Here we present the results of smoothed particle hydrodynamics simulations of high- 
velocity (up to ~ 5u esc ) collisions between 1 and 10 planets of initially terrestrial composition to 
investigate the end stages of formation of extrasolar super-Earths. As found in previous simulations 
of collisions between smaller bodies, when collision energies exceed simple merging, giant impacts are 
divided into two regimes: (1) disruption and (2) hit-and-run (a grazing inelastic collision and projectile 
escape). Disruption occurs when the impact parameter is near zero, when the projectile mass is small 
compared to the target, or at extremely high velocities. In the disruption regime, we derive the criteria 
for catastrophic disruption (when half the total colliding mass is lost), the transition energy between 
accretion and erosion, and a scaling law for the change in bulk composition (iron-to-silicate ratio) 
resulting from collisional stripping of a mantle. 

Subject headings: planets and satellites: formation — planetary systems: formation 



1. INTRODUCTION 

To date, more than 300 extrasolar planets have been 
discovered. Of these, more than 1 have masses < 10 M m . 
With t he launch of the CO ROT (|Borde et al.ll2003[) and 
Kepler (|Borucki et al.|[2003T ) satellites, it is expected that 
many transiting super-Earths will be discovered in the 
next few years. In fact, the first transiting super-Earth 
candidate, CoRoT-7b, has been announced (Leger et 
al, in prep). The radius determined from the transit, 
when combined with radial velocity measurements of the 
planet's mass, can be used to dete rmine a mean densit y 
and thus infer a bulk composition ([Valencia et al.ll2007f ). 

Although the chemical composition of planets in the 
solar system generally reflects the gradient in the neb- 
ula, the final collisions forming each planet may involve 
embryos scattered from regions with different bulk com- 
position; hence the final composition of each planet is 
thoug ht to be domin a ted by the last f ew impact events 
(e.g., IWetherilll 11994 IChambersI 12004( 1. Planet forma- 
tion simulations have found that late stage impacts can 
reach velocities up to six times the mutual escape ve- 
locity hi_the_te£restnal planet region of the solar sys- 
tem (jAgnor et al.l 119991 ). Due to the presence of gi- 
ant planets and planetary migration, impact velocities 
may be even higher in other planetary sy stems (e.g., 
iRavmond et~aTJ l2005bl : fMandell et al.ll2007t l. Thus, the 
formation of super-Earths may involve very high velocity 
impacts between large bodies (e.g., 0.1 to several Earth 
masses). Collision outcomes will vary widely depending 
on the impact velocit y, mass ratio between the bodies, 
and the impact angle (I Agnor k Asphaugll2004t I Asphaud 
l2009HStewart k Leinhardtll2009h . The conditions for im- 
perfect merging are particularly interesting as they pro- 
vide an opportunity to further alter the bulk composi- 
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tions of the growing planets. An erosive impact event 
may be responsible for the high bulk density of Mer- 
cury rel ative to the other terres trial planets in the solar 
system (|Benz et al.lll988l [2007( l. In addition, collisional 
erosion has been consid ered for the Earth-Moon impact 
(jO'Neill k Palmd[20ol) . 

The end stages of planet formation are typically stud- 
ied using iV-body codes that assume per fect merging 
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of th e masses of colliding e mbryos (e.g., 
19991 [Chambers! Isubmittedl ) . However, 
(1999) showed that the assumption of perfect merg- 
ing often leads to the combined body rotating faster 
than the critical limit. Recently, various hybrid tech- 
niques combine TV-body dyn amics and a statistical treat- 
ment of fragmentation (e. g . , ICharnoz" k Morbid clli 2003; 
iBromlev k Kenvonll2006( ). The adopted fragmentation 
models are usually derived from simulations of ero- 
sive collisions over a range of target sizes and im- 
pact velocities (e.g.,|Benz k Asphaugj[l999t iBenzl l2000t 
IStewart k Leinhardt I2009D . Models of planetary accre- 
tion (still assuming perfect merging) now encompass the 
formation of extra solar systems with sup er-Earths and 
giant planets (e.g.. iThommes et al.l l2008). Hence, there 
is a need for simple descriptions of the outcome of colli- 
sions between large planetary embryos. 

To complicate matters, recent studies of collisions be- 
tween planetary embryos describe two dramatically dif- 
ferent re gimes when the impact energies exceed perfect 
merging (I Agnor k Asphaud l200l lAsphaug et alJ 120061 : 
I Asphaud 1200911 : (1) erosion via disruption and gravita- 
tional dispersal and (2) "hit-and-run" events where the 
projectile grazes the target, which is left nearly intact. 
The boundary between the two regimes depen ds on the 
mass ratio and impact angle (see Figure 17 in lAsphaud 
2009]). The disruption regime is typically characterized 
by the catastrophic disruption criteria, Q* , which is a 
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size-dependent curve describing the projectile kinetic en- 
ergy divided by the target mass such that half the tar- 
get mass remains. For small bodies (< 1 km in diam- 
eter), tensile strength, which decreases with size, deter- 
mines the outcome. For larger bodies, gravity dominates 
and the largest remnan t includes reaccumulated material 
l|Melosh fc Rvanl Il997f ). The transition between accre- 
tion (net growth of the target) and erosion is determined 
by the dependence of t he mass of the largest re mnant 
on the impact energy. IBenz fc Asphauel ( 19991 ) simu- 
lated colliding rocky and icy bodies to determine the 
catastrophic disruption criteria from scales of centime- 
ters to hundreds of kilometers. They found that the ra- 
tio of the largest remnant mass, Mi r , to the target mass, 
Mtarg, depends linearly on the impact energy divided 
by target mass, Q, scaled to Q*. This simple relation- 
ship has been confounded by the fact that Q* depends 
on the target size, impact velocity, mass ratio, and ma- 
terial properties. Thus, previous work relied upon the 
availability of catastrophic disruption criteria for specific 
impact scenarios. However, Ste wart fc Leinhardtl (|2009l ) 
have recently developed a more general expression for the 
disruption criteria. The ratio Mi r /M tar g is replaced by 
Mi r /Mtot, where M to t is the total colliding mass. The 
impact energy divided by target mass, Q, is replaced by 
the reduced mass kinetic energy scaled to the total col- 
liding mass, Qr — h/jy? /M tot , where /x is the reduced 
mass and Vi is the impact velocity. Using these variables, 
the new catastrophic disruption criteria, Q* RD , at which 
Mi r /Mtot = 0-5 is not dependen t on th e mass ratio. Ad- 
ditionally, Stewart fc Leinhardtl (|2009f ) provide parame- 
ters that describe a wide range of material properties and 
impact velocities. However, to date, the catastrophic dis- 
ruption criteria has not been derived for bodies larger 
than 100 km. 

Here, we simulate a wide range of impacts between 
planetary embryos up to and including super-Earths to 
derive the conditions for accretion versus erosion. We 
consider bodies with initial compositions similar to Earth 
and derive scaling laws for catastrophic disruption and 
changes to the iron-to-silicate mass ratio resulting from 



erosion. 



2. METHOD 



2.1. Smoothed Particle Hydrodynamics (SPH) 

Hydrodynamic codes have been used extensively to 
simulate giant impacts. SPH is a Lagrangian technique 
in which the mass distribution is represented by a finite 
set of particles ev olved with time (jGingold fc M onaghan 
ll977tlLucvlll977f) . To date, SPH has been us ed to study 
giant impacts in th e early histo ries of Earth (|Benz et alj 
119861: ICanud l200l. M ercury (IBenz et all [19881. 1 20071) . 
Mars (|Marinova et al.ll2008f) . and Pluto (jCanurj [200^ . 
These studi es have all utilized an SPH code descended 
from that of IBenz et all (|1986| ). 

Here we present the results of simulati ons of giant im - 
pacts using the SPH code GADGET ([Springe!! 12001 . 
which has been tested and utilized extensively in cos- 
mology and galactic dynamics. In particular, GAD- 
GET employs a tree- based scheme for comp uting the 
self-gravity of objects (|Hernquist fc Kat9ll989l ). and, un- 
like other implementations of SPH, a fully conservative 
approach for integrating th e hydrodynamic equations 
(Springel & Hernquist 2002) that maintains energy and 



entropy co nservation even when smoothing leng ths vary 
adaptivelv (|Hernquistl[l99l ; lO'Shea et al.ll2005f) . 

We have modified GADGET to compute thermody- 
namic quantities by interpolating between elements of 
tabulated equations of state. The tables were gener- 
ated using a revised versio n of the semianalytic equa- 
tion of state m odel ANEOS ([Thompson fc Lausonlll972l : 
iMelosbl [2007L We used ANEOS parameters for Si0 2 
(jMeloshl 120071) . forsterite, and iron (Melosh, private 
communication). The EOS tables were finely gridded 
to resolve phase bou ndaries (in a manner similar to 
ISenft fc Stewartl l2008l appendix). The modified code 
was tested with three-dimensional simulations of impacts 
between plates of identical material at multiple resolu- 
tions and shock pressures. The calculated peak pres- 
sures agree w ith the analytic impedance match solution 
( Melosh 1989, p. 54). We performed two additional tests 
of the code, which are describ ed below. As in p revious 
studies of giant impacts (e.g., IBenz et ail 119861 ). shear 
and tensile strength are neglected because gravity is ex- 
pected to dominate in this size regime. 

2.2. Test One: Shock Pressure Decay 

When a body impacts a half space, beyond an initial 
isobaric region near the impact point, the peak shock 
pressure decays exponentially with distance, d, via 

(1) 
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where n is a function of impact velocity. Using a two- 
dimensional Eulerian hydrocode and a range of materials, 
iPierazzo et "all (|1997l ) found that n = (-1.84 ± 0.17) + 
(2.61 ± 0.14) log(Ui), where V. is the impact velocity in 
km s -1 . 

Using our modified v ersion of GADGET, we repro- 
duced the simulations of IPierazzo et all (1997) for Si02 
and iron. We modeled the pressure decay profile from 
a 10 km diameter impactor striking a large flat surface 
(essentially a half space) at speeds of 10, 20, 30, 40, and 
60 km s" 1 , shown in Figure [Ij,. We fit the exponent n 
with (-1.70 ± 0.33) + (2.45 ± 0.24)log(U l ), in excellent 
agreement with previous work (Figure [TJd) . 

2.3. Test Two: Moon-Forming Impact 

The iron depletion in the Moon, isotopic sim- 
ilarities between Earth and Moon, and the large 
angular momentum of the Earth-Moon system led 
to the gi ant impact hypothesis f or the origin of 
the M oon (jHartmann fc Davisl 119751 : [C amcron fc Ward 
119761) . SPH has been extensively employed in model- 
ing potential Luna r- forming impacts, most recently by 
ICanud (l200l[200l) . 

We repeate d two of the simulations from Table 1 of 
ICanud (|2004l ). We started from the same initial condi- 
tions: a target and impactor of 70% silicates (forsterite) 
and 30% iron, hydrostatic pressure profile, and an initial 
isentropic temperature profile with a surface tempera- 
ture of 2000 K. The target and impactor were allowed to 
settle to negligible particle velocities in isolation. The fi- 
nal masses of the proto-Earth and proto-Lunar disk were 
calculated using the iterative p rocedure described in the 
appendix of lCanup et al.l (|2001f ). Tabled] presents our re- 
sults from these two si mulation s, which are in excellent 
agreement with those of ICanupl i 
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2.4. Present Work 

We performed a series of more than 60 simulations of 
impacts onto super-Earth sized bodies. The initial con- 
ditions were targets and impactors of roughly terrestrial 
composition, comprised of 67% forsterite and 33% iron. 
The bodies were initialized with hydrostatic pressure pro- 
hies and temperature profiles calculated f or super-Earths 
with surface temperatures of about 350 K (Valencia ct al. 
l2006t Figure 4). Previous work di d not show a stro ng de- 
pendence on initial temperature (|Benz et al.| [2007). The 
bodies were allowed to settle to negligible particle veloc- 
ities in isolation with the temperature fixed. We used 
three different target masses, 5M®, and WMq, 

with impactors of 1/4, 1/2, and 3/4 the target mass. All 
the initial experiments employed an impact parameter 
of zero, and a subset was repeated with impact angles 
of 45°. A few cases of undifferentiated, pure forsterite 
bodies were considered as a more direct comparison to 
previous work on catastrophic disruption, which has not 
included differentiated bodies. The number of particles 
ranged from 1.25 x 10 5 to 7 x 10 5 (30-45 particles per 
target diameter), and the results were checked for sensi- 
tivity to resolution. 

The mass of the largest remnant following each impact 
was determined using a f ragment search algorithm sim- 
ilar to that employed by iBenz fc Asphaual ([1999) . The 
potential and kinetic energies of all particles were calcu- 
lated with respect to the particle closest to the potential 
minimum. The center of mass position and velocity were 
then calculated for all bound particles and the process 
was repeated on the remaining unbound particles until 
the calculation converged (typically after four or five it- 
erations). 

3. RESULTS 

Our results fall into two regimes: (1) disruptive events 
and (2) hit-and-run events. At impact energies just 
above merging for events with projectile to ta rget mass 
ratios ranging from 1/10 to 1, the results from I Asphaud 
( 2009) indicate that the onset of hit-and-run events oc- 
curs at impact angles between 30° and 45° (where 0° is 
head-on). Our results are in agreement and all the 45° 
impact cases were merging or hit-and-run events. 

3.1. Catastrophic Disruption: Scaling largest remnant 
and iron composition 

In the disruption regime, the outcome of a collision is 
described in terms of the largest remnant mass and the 
catastrophic disruption criteria. In Figure [21 we present 
the largest remnant mass as a function of the scaled im- 
pact energy (diamonds) for our head-on impact events 
(open symbols) and a subset of sim ulations from can- 
didat e Mercury forming impacts from lBenz et al.1 (|1988l 
I2007h (filled symbols). The results are in excellent agree- 
ment with th e universal law for the l argest remnant mass 
presented by I Stewart fc Leinhardtl (|2009l ) (dotted line), 
where 

M lr /M tot = -0.5(Qr/Q*rd - 1) + 0.5 . (2) 

This l inear relationship derived bv I Stewart fc Leinhardtl 
(2009) for smaller bodies and slower impact velocities 
continues to hold into the super-Earth regime. Note 
that the utility of equation [2] is bounded by Qr/Q* rd 



of (events near merging) and 2 (a super-catastrophic 
regime). 

We find that the iron mass fraction of the largest rem- 
nant follows a simple power law (circles in Figure [2]), 
given by 



MFe/M w = 0.33 + 0.25{Q R /Q* RD ) 



1.65 



(3) 



Note that we included only the subset of simulations 
from Benz from which the catastrophic disruption cri- 
teria could be derived. 

The fitted catastrophic disru ption criteria from our 
simulations and the subset from lBenz et a l. (1988. f2007h 
are shown in Figure [3] The lines are the gravity-regime 
velocity-dependent cata strophic disruption c r iteria for 
low strength bodies from lStewart fc Leinhardtl (|2009f ) (in 
cgs units): 

Q*rd = io- 4 RhiK 08 - (4) 

Rci is the radius of a spherical body with M to t and den- 
sity of 1 g cm~ 3 . The Q* RD data points for giant impacts 
fall short of the prediction from Equation [4] by less than 
a factor of 2. Hence, the results for super-Earths are 
entirely consistent with disruption criteria developed for 
smaller bodies in the gravity regime, and one can rea- 
sonably determine the largest remnant mass and bulk 
composition following a collision knowing only the total 
colliding mass and impact velocity. 

3.2. Disruption Versus hit-and-run impacts 

Our high velocity impacts at angles of 45° by projectile 
masses between 1 /4 and 3/4 the target mass produced no 
significant change to the mass or the bulk composition 
of the target. In such cases, there was only an inelas- 
tic collision resulting in the target heating and acquiring 
angular momentum. Hence, we independ ently verify the 
pheno menon described as hit-and-run bv lAsphaug et aTl 

pool . 

The contrast between the disruption and hit-and-run 
re gimes is il l ustrat ed by the accretion efficiency defined 
bv I Asphaud (120091 ) . 



M, r - Mi 



targ 



Al, 



(5) 



proj 



where positive values correspond to growth and negative 
values to erosion of the target. Note that the value of £ 
that corresponds to catastrophic disruption depends on 
the projectile to target mass ratio (— 0.5M ta rg/M pro j + 
0.5). Figure |H presents accretion efficiency as a function 
of impact velocity, scaled to the mutual escape velocity 
{v 2 esc = 2G{M targ + M proj )/(R targ + Rproj))- For head- 
on impacts (filled black symbols), the fraction of the col- 
liding mass remaining in the largest remnant decreases 
steadily with increasing Vi/v esc . However, for 45° im- 
pacts (open symbols), the outcome is nearly a step func- 
tion from perfect merging to remnant masses approxi- 
mately equal to the target mass (£ = 0). Our results for 
super-Earths are in excellent agreement with the studies 
of accretion efficiency for i mpacts onto much smaller bod- 
ies (O.lOMm) pr esented in lAgnor fc Asphaud (120041 ) and 
lAsphaud (120091 ) (e.g., dotted line in Figure Ha). Note 
that the hit-and-run regime encompasses a very large 
range of impact velocities. In this work, we considered 
events exceeding 5i> esc , about twice as fast as in previous 
work. Hence, the possible outcomes of impacts between 
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comparably sized bodies are similar over orders of mag- 
nitude in target mass. 

At sufficiently high impact energies, the hit-and-run 
regime transitions to disruption. A study of collisions 
onto 100 km diameter bodies by projectile masses up to 
a ratio of 1:10 found that the impact energies required 
for catastrophic disruption may increase by an order of 
magn itude or more wit h increasing impact angle (Figure 
1 in iDurda et a l. 2007). The impact angle dependence 
is seen in Figure 2)d: near the catastrophic disruption 
value of £ = —2 for 1:5 mass ratio (stars), compare the 
offset in velocity for impact angles of 0°, 30°, and 45°. 
The increased disruption energy as a function of impact 
angle is given in the diamond points in Figure [31 

Thus, the results of giant impact events at velocities up 
to several times v esc depend primarily on (1) the impact 
angle and mass ratio which define the transition between 
the disruption and hit-and-run regimes and (2) the total 
colliding mass and impact velocity which define Q* RD in 
the disruption regime. The outcome falls in the disrup- 
tion regime at any impact angle when the projectile is 
very small and fast but only when the impact angle is 
small (typically < 30°) when the projectile and target 
are comparably sized. Because 45° is the most probable 
impact angle, catastrophic disruption of protoplanets is 
most likely to occur when the impactor is much smaller 
than the target. 

4. DISCUSSION AND CONCLUSIONS 

Using the results from this work and previous studies 
of giant impacts, we have described the outcome of col- 
lisions between a wide range of masses, from planetary 
embryos to super-Earths . In the disruptive regime, the 
framework developed by I Stewart fc Leinhardtl (|2009[ ) to 
describe the transition between accretion versus erosion 
in hypervelocity impacts holds from the lab scale, with 
sizes in centimeters and velocities in m s -1 , up to the 
scales of the largest expected super-Earths, with sizes up 
to ten thousand kilometers. 

In the disruption regime, we derived a scaling law for 
the changes in bulk composition resulting from collisions 
between differentiated bodies of Earth-like composition. 
Based on the boundary between disruptive and hit-and- 
run type events, significant compositional changes re- 
quire either a small impact angle or a small, fast projec- 
tile. Because of this narrow range of impact conditions 
that can lead to an increase in bulk density, the presence 
of unusually dense planets can place constraints on the 
dynamical history of a planetary system. In particular, 
if planet growth occurs by merging similar-sized bodies, 
changes to the bulk composition of the target are un- 
likely. Thus, one-time high-energy impacts provide the 
best means of increasing a planet's bulk density. If the 
specifics of the planetary system rule out such a colli- 
sion, then such a planet is more likely to have formed in 
a protoplanetary disk with substantially different chem- 
ical abundances or with a density profil e quite different 
from t he minimum-mass solar nebula (jRavmond et al.l 
l2005aft . 

For very energetic impacts, the largest remnant can 
have a substantially different composition than that of 
Earth, even when the target and impactor initially had 
Earth- like bulk compositions. This is consistent with the 
hypothesis that Mercury, which is ~ 70% iron by mass, 



owes its present bulk composition to a catastrophic im- 
pact event. For planet masses close to that of Mercury, 
the critical velocity needed to yield a remnant with such 
an iron enrichment is abou t 10-20 km s _1 , as was found 
bv iBenz et alJ (|1988l 120071 ). Such velocities are smaller 
than the orbital velocities inside 1 AU for a ~ 1M Q star. 
However, for more massive planets, the critical veloci- 
ties can enter the 50-60 km s -1 range, requiring extreme 
conditions for a Mercury-like planet to be created. Thus, 
Mercury-like planets should be much more common than 
super-Mercuries, although even these may form in very 
close orbits in extrasolar systems. Hence, the probabil- 
ity of forming super-Mercuries warrants further investi- 
gation (Marcus et al., in prep). 

We thank Joe Barranco for help in code development 
and Zoe Leinhardt, T.J. Cox and the anonymous re- 
viewer for useful feedback. The simulations in this Letter 
were run on the Harvard FAS Odyssey cluster. 
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TABLE 1 

Test Two: Potential Moon forming impacts. 



N/W 3 7 v b M D /M L L D /L EM M Fe /M D M imp /M D L F /L EM M M /M L 



This work 



120 


0.13 


1.00 


0.730 


1.40 


0.27 0.04 


0.78 


1.06 


1.24 


120 


0.13 


1.00 


0.740 


1.44 


0.30 0.04 


0.78 


1.07 


1.50 












Canup (2004) 








120 


0.13 


1.00 


0.730 


1.59 


0.30 0.05 


0.77 


1.16 


1.28 


120 


0.13 


1.00 


0.740 


1.54 


0.31 0.04 


0.80 


1.19 


1.44 



Note. — N — number of particles, 7 — mass ratio, v - impact velocity in units of mutual escape 
velocity, b — impact parameter, Mp, - disk mass, - Lunar mass, Lp> - disk angular momentum, L E m 
- angular momentum of Earth-Moon system, Mp e - mass of iron in disk, Mi mp - mass of disk originating 
from impactor, hp - final system angular momentum, and Mm - estimated mass of satellite forming from 
disk. 
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Fig. 1. — Test One: Shock Pressure Decay, (a) Shock pres- 
sure as a function of scaled distance from impact point for a SiC>2 
target and projectile, (b) Power-law decay exponent (equation [TJ 
as a function of imp act velocity, with fits from present work and 
IPierazzo et~aTI (1997ft . 




Fig. 2. — Largest remnant mass and iron fractio n vs. scaled 
impact energ y. Open symbols, this work; filled symbols. Be nz et al.l 
il988l 120071) . Dotted line, Equation dashed curve, Equation [3] 




Fig. 3. — Catastrophic disruption criteria Q^jj vs. RCi. Re- 
sults from impacts between differentiated bodies (x) and pure 
forste rite bodies (squares). Diamonds are Mercury-size events 
from IBenz et al.l p988l . |200W) where shading indicates 0°, 30° 
and 45° impacts. Symbol colors indicate impact velocity. Ver- 
tical black bars are la errors on Q^jj- Solid lines are disrup- 
tion criteria for impact velocities of 10 (black), 20 (blue), and 40 
km s -1 (orange) (Equation |4). The inset shows full size range 
for the velocity-dependent catastrophic disruption curves from 
IStewart fc Leinhardtl l(2009h . 
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Fig. 4. — Accretion efficiency, £, as a function of scaled im- 
pact velocity, (a) Accretion regime. Collision outcome depends 
strongly on impact angle . Dashed l ine pr esents results for 45°, 1:2 
mass ratio impacts from [Asphaug (2009, Figure 17). (b) Disrup- 
tion regime. Catastrophic disruption is § = —4.5, -2.0, and -0.5 for 
1:10, 1:5, and 1:2 mass ratios, respectively. In both panels, sym- 
bols denote mass ratio and shading denotes imp act angle. Circles, 
pentagons, stars - from Bcnz et al. (1988, 2007). Other symbols - 
this work. 
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